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We present a Full Potential Multiple Scattering (FP-MS) scheme for the interpretation of several 
X-ray spectroscopies that is a straightforward generalization of the more conventional Muffin- Tin 
(MT) version. Like this latter, it preserves the intuitive description of the physical process under 
consideration and overcomes some of the limitations of the existing FP-MS codes. It hinges on a 
fast and efficient method for solving the single cell scattering problem that avoids the convergence 
drawbacks of the angular momentum (AM) expansion of the cell shape function and relies on an 
alternative derivation of the multiple scattering equations (MSE) that allows us to work reliably 
with only one truncation parameter, i.e. the number of local basis functions in the expansion of the 
global scattering function determined by the classical relation Z max ~ k R. 

PACS numbers: 61.05.js,61.05.cj,71.15.Ap 
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Multiple scattering theory ( MST ) has been widely 
used to solve the Schrodinger equation (SE) ( or the as- 
sociated Lippmann-Schwinger equation (LSE)) both for 
scattering and bound states. It was proposed originally 
by Korringa and by Kohn and Rostoker (KKR) as a con- 
venient way to calculate the electronic structure of solids 
[H, 0] with potentials of the muffin-tin (MT) form (i.e. 
potentials that arc bounded by non overlapping spheres 
and spherically symmetric), and was later extended to 
polyatomic molecules by Slater and Johnson How- 
ever the MT approximation cannot properly describe a 
great number of physical systems, ranging from open lat- 
tices to molecular systems with substantial anisotropy ( 
e.g. systems of biological interest ), to surfaces and in- 
terfaces. 

The attempts to extend MST to space-filling potential 
cells in order to eliminate the interstitial region and take 
full account of the asphericity of the potential, have gen- 
erated a lot of controversies that have gone on for more 
than twenty years Many (but not all) of the ques- 
tions are now settled and we refer the reader to the book 
of Gonis and Butler [|| for a comprehensive review of 
the state of the art in this field. Usually the applica- 
tions of the space-filling method have regarded mainly 
the calculations of the electronic structure of solids, i.e. 
states below the Fermi level. Applications to states well 
above the Fermi energy, as required in the simulations 
of x-ray spectroscopies, like absorption, photo-emission, 
anomalous scattering, etc., have been scarce, probably 
due to the difficulties encountered in the numerical im- 
plementation of the method. We mention here the work 
by Huhne and Ebert [f| on the calculation of x-ray ab- 
sorption spectra using the full-potential spin-polarized 
relativistic MST that of Ankudinov and Rehr [?J and 
that of Foulis [8| based on a version of the MST that 
uses spherical cells and treats the interstitial potential in 
the Born approximation 0, ■ 

All these methods, however, have their limitations and 



drawbacks. The method used by Foulis treats in an ap- 
proximate way the potential in the interstitial region and 
moreover looses one of the major advantages of the MST, 
namely the separation between dynamics and geometry 
in the solution of the scattering problem. Huhne and 
Ankudinov use the potential shape function to gener- 
ate the local basis functions which are at the heart of 
MST. The expansion of the shape function and the cell 
potential in spherical harmonics leads to a high number 
of spherical components in the coupled radial equations 
that becomes progressively cumbersome to handle and 
time consuming with increasing energy and in absence of 
symmetry. This feature might also be at the origin of an- 
other problem related with the saturation of "internal" 
sums in the MSE Q, as discussed below. Moreover no 
critical discussion is devoted in their work to the conver- 
gence problems of MST. 

The purpose of this letter is to give an alternative 
derivation and interpretation of the FP-MS equations 
that will allow us to work with square matrices for the 
phase functions Sll' and Et,l' and for the cell 111' 
matrix (see below) with only one truncation parameter, 
contrary to the present accepted view [f|. As a result 
this scheme can be viewed as a natural extension of its 
MT counterpart, with all the consequent advantages for 
the interpretation of X-ray spectroscopies. In connection 
with this we shall present a new scheme to generate local 
basis functions for the truncated potential cells that is 
simple, fast, efficient, valid for any shape of the cell and 
reduces to the minimum the number of spherical harmon- 
ics in the expansion of the scattering wave function. 

In order to solve for scattering states we seek a so- 
lution of the SE continuous in the whole space with its 
first derivatives, satisfying the asymptotic boundary con- 



dition V(r;k) ~ (^) 



/(*;k)' 



where k 



is the photo-electron wave- vector. We partition the space 
in terms of non overlapping space-filling cells f2j with 
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surfaces Sj and centers at R 3 . This is equivalent to par- 
titioning the overall space potential V(r) into cell poten- 
tials, such that V(r) = J2j v j( r j), where Vj(rj) takes the 
value of V(r) for r inside cell j and vanishes elsewhere. 
Here and in the following rj = r Rj. The partition 
is assumed to satisfy the requirement that the shortest 
inter-cell vector R,j = R; — Rj joining the origins of the 
nearest neighbors cells i and j, is larger than any intra- 
cell vector or rj, when r is inside either cell i or cell 
j. We also assume that there exists a finite neighbor- 
hood around the origin of each cell lying in the domain 
of the cell Q . We then start from the following identity 
involving surface integrals in dr = da 
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j [G+(r' - r)VV(r) - </>(r)VG+(r' - r)] • n, d^ 

= j [G+ (r' - r) W(r) - ^(r) VG+ (r' - r)] • n da . 

Here Gq(t' — r) is the free Green's function and £! = 
■ f2j, with surface S . This identity is valid for all r' 
lying in the neighborhood of the origin of each cell, since 
in this case the integrands are continuous with their first 
derivatives. 

The heart of MST is the introduction of the func- 
tions $i(rj) which inside cell j are local solutions of 
the SE with potential Vj(rj) behaving as Ji(rj) for 
rj — ► 0. They form a complete set of basis functions 
such that the global scattering wave function can be lo- 
cally expanded as tp{rj) = ^^W^ifc) El- Us- 
ing this expansion in the above surface integrals, taking 
r' in the neighborhood of the origin of cell i, so that 
G+(r' - r) = G+(r; - r<) = El MO^lfr) ( since r » 
confined to lie on the cell surfaces), and putting to zero 
the coefficients of Jl(i^) due to their linear independence, 
we readily arrive at the MST compatibility equations for 
the amplitudes ^(k) (see Ref. [B[ page 129 for analogous 
derivation in the case of bound states) 

Y,H l i L ,Al,(k) =Y L (ky k - R '(k/n)^ = F L (k) (1) 
where 

H% = [ [H+ (t<) V$ L , fa ) - $ L , (tj ) VH t (r, ) ] ■ n 3 da, 

and the rhs term comes from the outer sphere Q (see 
Appendix A of Ref. [ll[ for details). As usual Ji,(r) = 
ji(kr)Y L (r) and H+(r) = -ikhf (kr)Y L (r) . The usual 
derivation of the MSE now proceeds by re-expanding 
H~i(vi) around center j by use of the equation iJ^(r.j) = 
El 1 G*l L iJL<(rj), where G % [ v are the free electron prop- 
agator in the site and angular momentum basis ( KKR 
real space structure factors) . Unfortunately this relation 
introduces a further expansion parameter into the theory 



(with related convergence problems) which is actually un- 
necessary, as shown below. 

We in fact observe that the integrals over the surfaces 
of the various cells j can be calculated over the surfaces of 
the corresponding bounding spheres (with radius R ] b ) by 
application of the Green's theorem, since both -ff^ (r) and 
$L(r) satisfy the Helmholtz equation outside the domain 
of the cell (where the potential is zero) . We then use the 
following relation 



F^(f,)V7J+(r 4 ) -njdff,- = Gf L , —jjv(kRi) (2) 

dR l 



(and the similar one without derivatives) which is 
exact for all L provided |r, — Tj\ = Rij > rj 
for r lying on the surface Sj. This is a conse- 
quence of the fact that under this condition the series 
Hi(ri) = El' ^LL"^ L '( r i) converges absolutely, since 
G l l L ,Juivj) < {rj/Rijf'l/ikRij^W + + for 
fixed I, as can be seen by using the usual expression for 
G l [ L , and the asymptotic expansion of the Bessel func- 
tions for high values of the index V (V ^> krj) [12]. By 
use of the Weierstrass critcrium, the series is also uni- 
formly convergent in the entire solid angle domain and 
can therefore be integrated term by term (l3j , leading to 
the desired result (this property is also true for the series 
derived with respect to r) . 

By inserting in Eq. ([1]) the expression for the ba- 
sis functions expanded in spherical harmonics $z,(r) = 
J2l< Rl'l{i~)Yl' (f) and using Eq. ©, we finally obtain 



Y,E i LL/ A L ,(k) - J2 G i i v ,si, lv Ai,(k) = ii(k) 

L> j,L',L" 

(3) 

defining 

E LL > = R 2 b W[-ikh+,R LL ,}; S LU = R 2 b W[j u R LL ,] 

where the Wronskians W[f, g] = fg' — gf are calculated 
at Rt and reduce to diagonal matrices for MT potentials. 
We observe that, even if the potential has a step, the 
wave function and its first derivatives are continuous, so 



that the AM expansion converges uniformly in f 14| and 
can be thus integrated term by term. Eq. ([3]) looks for- 
mally similar to the usual MSE. However we notice that 
the sum over L" runs over the angular momentum com- 
ponents of the basis functions $ l (r) and is not affected 
by convergence constraints related to the re-expansion of 
Hi(r.i) around center j. 

To find the local solutions of the SE we do not expand 
the truncated cell potential to avoid AM expansion prob- 
lems. Dropping the cell index j, we write the SE in polar 
coordinates for the function Pl(?) = r^^i^) 



d 2 . M 

— + E-v(r,r) 



P L {r 1 v) = ^L 2 P L {r 1 v) (4) 
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where L 2 is the angular momentum operator, whose ac- 
tion on Pt(r, f) can be calculated as: 

L 2 P L {r, r) *'(*' + l)rR L <L(r)Y L ,(r). (5) 



Equation in the variable r looks like a second order 
equation with an inhomogeneous term. Accordingly we 
use Numerov's method to solve it. As is well known, 
putting = Pj,(rj,fj-) and dropping for simplicity the 
index L, the associated three point recursion relation is 

where, 



^*>i — 1 Y2 Vi '^ 



Bi 



5h/ 
~6 



vn = 12 - 10A, 



= v(n,Tj)-E 
h 2 

9i,j = ^2 + 10 9i,j + Si-ij] 



Here i is an index of radial mesh and j an index of an- 
gular points on a Lebedev surface grid [l5j|- Obviously 
RL'h{ri) = Y J j w o p L{ri,rj)/r i Y L >{rj), where uij is the 
weight function for angular integration associated with 
the chosen grid. 

Only the inhomogeneous term qi+ij in the recurrence 
relation, containing the still unknown term fi+i.j, pre- 
vents us to solve the equation by iteration, from the 
knowledge of fij and fi-i.j at all the angular points. 
This difficulty is easily overcome by introducing the back- 
ward second derivative formula, whereby 
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h 5 

[lSqtj - 2ft_xj + qi-2,j] + 
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Qij (6) 



so that at the cost of a small errors 0(h 5 ) only the back- 
ward points fi_j, fi—i,j and are now involved. The 
appearance of q"', strictly infinite at the step point, does 
not cause practical problems. 

In this way the three-dimensional discretized equation 
can be solved along the radial direction for all angles in 
an onion-like way [12| , provided the expansion flS)) is per- 
formed at each new radial mesh point. We use a log-linear 
mesh p = ar + /31nr, to reduce numerical errors around 
the origin and the bounding sphere. [l6|. We tested this 
modified Numerov method against analytically solvable 
separable model potentials, with and without shape trun- 
cation, obtaining very good results (l2j . 

The next delicate point to tackle is how to handle and 
truncate the various L sums in Eq. ([3]). Here only two 



truncation parameters appear, the number of basis func- 
tions and the number of their AM components, corre- 
sponding to the indexes V and L in Rll 1 - (The exter- 
nal index L in Eq. |T|). coming from the expansion of 
Gq, must coincide with the index L of Rll' when cal- 
culating H" L , = E % LL ,, due to the orthogonality of the 
spherical harmonics). These two indexes are in principle 
unrelated, although one can speculate that for I > k Rb 
and positive energies the wave-function hardly sees the 
anisotropy of the potential. This is in practice what hap- 
pens in our method of generating the local basis func- 
tions, so that the two indexes can be cut safely at the 
same value i max ~ k Rb and the matrices S and E can 
be treated like square matrices. This is also in keeping 
with the physical fact that the elements of the atomic T- 
matrix Tu = (SE~ 1 )u tend rapidly to zero for I > l max . 
The need to converge first the internal sum over L" in Eq. 
([3]), pointed out by various authors, was probably related 
to the slow rate of convergence in the AM expansion of 
the basis functions, due to their method of generation. 

We are aware that even in the case of "small" over- 
lap of the cell bounding spheres the L truncation pro 
ccdurc is likely to be divergent 
tion 
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A simple calcula- 
V2\ shows that G{R lj ) a T u goes like {2R h / R^) 21 / I s 



for I 



ij III J-U B uco iJJVO V*-* 1 */ ) 

which explains why for "small" overlap 



the spectrum seems to converge at first, but actually di- 
verges for high values of l max , which in practice depends 
on various parameters, with the typical behavior of the 
asymptotic series. In keeping with this view one can 
attribute a meaning to the finite result obtained by trun- 
cating the MSE: at a given i max the difference between 
the exact and the approximate result can be made very 
small, provided the independent parameters of the theory 
(like the photo-electron energy and the amount of overlap 
between the bounding spheres) lie in a definite domain of 
their definition space. For energies in the near edge re- 
gion this can be obtained by taking a moderate overlap 
between bounding spheres, of the order of 30~40%, so as 
to reduce the space between them and their MT spheres. 
Empty spheres can be added to satisfy this condition. 

We have tested the present FP-MS scheme against 
the analytical solution of the absorption cross section for 
hydrogen-like atoms given by (l8j 



Cr(fc) = 47T 2 Q' 



2 _J_ 

3 Z 2 



-4f tan-^i) 



1 + (|) 2 



1 



-27T ^ 



in the case of the Li 2+ atom (Z = 3). Even though 
the potential is spherically symmetric in the whole space 
with respect to the atomic center so that it is easy to 
reproduce numerically the cross section, this is not ob- 
vious in the MS scheme. To this purpose we have par- 
titioned the space inside a sphere of radius R = 8.6 au 
into an atomic sphere of 4.15 au and 14 other empty 
spheres, all truncated so that the resulting polyhcdra do 
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FIG. 1: (color online) Cross section for Li 2+ with 15 cells 
compared to the analytical result. The solution for a MT 
central sphere is also shown. 
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not overlap and such that their bounding sphere do not 
overlap more than 40%. To calculate the contribution 
of the outer sphere we integrated inwardly the Coulomb 
potential. Fig. [T] shows the almost exact agreement be- 
tween the analytical and the numerical result, indicating 
that the partitioning procedure for solving the SE is able 
to reconstruct the global solution. Moreover the oscilla- 
tions due to the truncation of the potential inside each 
cell (shown by the solution for a truncated central sphere 
with radius 4.15 au) cancel each other, showing that at a 
common boundary the overall solutions inside two adja- 
cent cells are continuously smooth. For this test a value 
of Imax = 4.15 -\/3 ~ 8 was taken at the end of the energy 
interval E max = 3 Ryd. 

Fig. [2k shows an application of the method to the 
calculation of the Ge K-edge absorption spectrum of the 
tetrahedral molecule GeCLi [3]. The MT approxima- 
tion could never reproduce the first bump after the main 
transition. Its appearance is due to the introduction of 
the anisotropy of the potential inside the atoms and the 
presence of four empty Voronoi cells completing the BCC 
unit cell. An i max = 4 was sufficient to reach convergence 
of the spectrum, as verified by using higher I values up 
to l max = 10 (Fig. [lb). 

In conclusions we have developed a FP-MS scheme 
which is a straightforward generalization of the usual 
theory with MT potentials and implemented the code 
to calculate the cross section for several spectroscopies 
like absorption, photo-electron diffraction and anomalous 
scattering. The key point in this approach is the gener- 
ation of the cell solutions <f>£ ( r ) for a general truncated 
potential free of the well known convergence problems of 
AM expansion together with an alternative derivation of 
the MSE which allows us to treat the matrices S and E 
as square, with only one truncation parameter given by 
the classical relation Z max ~ k R. Even though this trun- 
cation procedure does not converge, taking a moderate 
overlap between bounding spheres assures satisfactory re- 
sult in the approximation of the exact solution. At the 
same time we have provided an efficient and fast method 



FIG. 2: (color online) (a) Cross section for GeCU molecule 
with nine scattering cells located at the sites of a BCC lattice, 
compared with the MT result and experiment, (b) Study of 
its convergence rate as a function of l ma x up to Imax ~ 10. 



for solving numerically a partial differential equation of 
the elliptic type in polar coordinates which can also be 
used to solve the Poisson equation. 
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